Disorder and interactions in quantum Hall ferromagnets near v — 1 
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We report on a finite-size Hartree-Fock study of the competition between disorder and interactions 
in a two-dimensional electron gas near Landau level filling factor v = 1. The ground state at v — 1 
evolves with increasing disorder from a fully spin-polarized ferromagnet with a charge gap, to a 
partially spin-polarized ferromagnetic Anderson insulator, to a quasi-metallic paramagnet at the 
critical point between i = and i = 2 quantum Hall plateaus. Away from v = 1, the ground state 
evolves from a ferromagnetic Skyrmion quasiparticle glass, to a conventional quasiparticle glass, and 
finally to a conventional Anderson insulator. We comment on signatures of these different regimes 
in low-temperature transport and NMR lineshape and peak position data. 
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I. INTRODUCTION 

Because of the macroscopic degeneracy of single- 
particle states in a Landau level, neither disorder nor 
electron-electron interactions in a two-dimensional elec- 
tron system (2DES) can be treated perturbatively in the 
quantum Hall regime. This is the raison dletre for the 
many interesting and surprising phenomenal!] which have 
arisen in quantum Hall physics. Theories of quantum 
Hall physics usually include either only interactions or 
only disorder, although both are always present. In par- 
ticular, it is common to include only disorder in studies 
of the integer quantum Hall effect (IQHE), which gen- 
erally focus on the sudden jump in the Hall conductiv- 
ity between values separated by e 2 //i, and it is common 
to include only interactions in studies of the fractional 
quantum Hall effect (FQHE), which generally focus on 
the ability of interactions to create charge gaps at partial 
Landau level fillings. The competition between interac- 
tions and disorder has often, but not always,B'cl been ne- 
glected, in part because of the lack of easily manageable 
analytical and numerical tools that can deal with both 
simultaneously. In this paper we address an instance in 
which this competition is particularly direct and can be 
successfully addressed with elementary techniques. 

At Landau level filling factor v = 1, the ground state 
of a disorder-free 2DES is a strong ferromagnet p i.e. it is 
completely spin-polarized by a Zeeman field of infinitesi- 
mal strength. In practice, of course, the field experienced 
by a 2DES in the quantum Hall regime is not infinites- 
imal; however the field's Zeeman coupling to the elec- 
tron spin is typically very weak compared to other energy 
scales. (In referring to these systems as ferromagnets we 
are emphasizing that they remain spin polarized in t 
limit of zero Zeeman splitting. In experimental systems,! 
typical values of the interaction and Zeeman energy scales 
are ~ 160 K and ~ 3 K respectively. The spin-splitting 
produced by this bare Zeeman coupling is usually negli- 



gible in paramagnetic states.) The quantum Hall ferro- 
magnet has a large gap for charge excitations, and hence 
has a robust quantum Hall effect. For typical Zeeman 
coupling strengths, its elementary charged excitations 
are topologically charged spin, textures (Skyrmions) con- 
taining several flipped spins .Bern Large Skyrmions have 
a Hartree energy cost smaller that those of conventional 
quasiparticles and, because the.-spms align locally, only 
slightly higher exchange energy.0 Because Skyrmions are 
the lowest-energy charged excitations, the global electron 
spin polarization is expected to decrease rapidly as |1 — v\ 
increases. This has been observed, in nuclear magnetic 
resonance (NMR) experiments-MB On the other hand, 
the v = 1 state of non-interacting disordered electrons 
differs qualitatively. The ground state is a compress- 
ible paramagnet with no Knight shift and no gap for 
charged excitations. For zero Zeeman coupling, quasipar- 
ticle states at the Fermi energy are quasi-extended and 
cause the Hall conductivity to suddenly jump by 2e 2 /h 
as this filling factor is crossed; v = 1 is in the middle 
of a Hall 'riser', instead of being at the middle of a Hall 
plateau. 

The competition between disorder and interactions at 
v = 1 can be addressed^ using the Hartree-Fock ap- 
proximation which has the virtue of being exactEJ in 
both the non-interacting and the non-disordered lim- 
its. Experimental information on this competition comes 
primarily from transport and NMR studies. Early 
NMR studies^ of weak-disorder quantum Hall ferro- 
magnets, yielded relatively featureless lincshapes and 
Knight shifts in SQftd, agreement with Hartree-Fock (HF) 
theory estimatesS'LrEl of Zeeman-coupling and filling- 
factor dependent Skyrmion sizes. (Effective field the- 
ory estimatesLrEHl are not accurate in the case of typical 
Zeeman coupling strengths.) More recent experimentsEJ 
paint a more complex picture, in part because the mea- 
surements were performed at lower temperatures where 
the signal is not motionally averag edHB It is now clear 
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that disorder plays a role in the interpretation of these 
experiments, even when it is weak. At stronger disor- 
der, as the non-interacting limit is approached, the spin- 
polarization must eventually vanish. In our model cal- 
culations we find that as the interaction strength is in- 
creased relative to disorder at v — 1, the 2DES ground 
state suffers a continuous phase transition a from para- 
magnetic to a ferromagnetic state. Depending on the 
details of the disorder model, a second continuous phase 
transition to a fully spin-polarized incompressible strong 
ferromagnet with a gap for charged excitations may oc- 
cur at still stronger interactions. For the disorder mod- 
els we use, the fully polarized state is reached when 
the Coulomb energy scale is approximately twice the 
Landau-level-broadening disorder energy scale. Away 
from v = 1, screening by mobile charges reduces the 
importance of disorder and the system reaches maximal 
spin-polarization at smaller interaction strengths. The 
maximally polarized ground state at moderate interac- 
tion strengths is best described as a glass of localized 
conventional quasiparticles formed in the v = 1 fully po- 
larized vacuum. Only for stronger interactions do we find 
a phase transition to a state with non-collinear magne- 
tization in which the localized particles have Skyrmionic 
character. 

We organize this paper as follows. In Sec. II we sum- 
marize our implementation of finite-size HF theory in 
the lowest Landau level (LLL). In Sec. Ill and IV we 
present and discuss our numerical results for calculations 
at v = 1 and v ^ 1. The possibility, discussed in re- 
cent work,B that at v = 1 disorder will induce reduced- 
size Skyrmion-anti-Skyrmion pairs in the ground state is 
specifically addressed in Sec. III. Finally, in Sec. V we 
present our conclusions. 



II. HARTREE-FOCK THEORY IN THE LLL 

The HF approximation allows the interplay between 
disorder and interactions to be addressed while retain- 
ing a simple independent-particle picture of the many- 
body ground state. For the current study, the use of 
this approximation is underpinned by the fact that it re- 
produces the exact ground state at v = 1 in both weak 
interaction and strong interactionEl limits. HF theory is 
a self-consistent mean field theory, and as such it has 
strengths and shortcomings, which we discuss later. In 
this section, we outline the basic formalism of HF ap- 
proximation calculations in the LLL limit. 

In a strong magnetic field, the Landau level splitting is 
very large and, since excitations to higher Landau levels 
are effectively forbidden at the low experimental tem- 
peratures, we follow the common practice of considering 
only LLL states. Neglecting the frozen kinetic-energy de- 
gree of freedom, the Hamiltonian in second quantization 
is written as 



where TLi is the interaction part of the Hamiltonian 
H!= l -JdvJ dr'^ % (r-r')^(r)^(r')^(r')^« 

a a' 

Tidis is the external disorder part of the Hamiltonian 
Hdis = / dr^u E (r)^+(r)Vv(r) , 

and Tiz is the Zeeman term 

■Hz = -\g^ B J dr^t,(rOVv(r')r CT v • B(r) , 

with a =\,\, vi and ve being the Coulomb interac- 
tion and disorder potentials respectively, and t\ being 
the Pauli matrices. Here we also define the Zeeman cou- 
pling strength as g — g^sB /(e 2 /el) for later reference. 
We chose the the Landau gauge elliptic theta functions 
as our basis 

where x m , s = + sL x , m, to' = 1, . . . , Nj,, N<p 

A (2?rZ 2 ), and I is the magnetic length which we set equal 
to 1 for simplicity. These wavefunctions satisfy the semi- 
periodic boundary conditions 4> m {x,y) — 4> m (x,y + L y ) 
and 4> m (x + L x ) = cx.p(+iL x y/l 2 )(j) m (x,y). 

We consider the HF single particle states to be a linear 
combination of the up and down spin states of these basis 
functions 

\a) = ^^(mer|a) \ma) . 

m,er 

Before writing down the Hamiltonian matrix in the HF 
approximation, we introduce several notation simplifying 
definitions, closely following previous HF studies B The 
expectation value of the particle density (in momentum 
space) is given by 

(p(q)) =Y,n F {e F - e a ){a\e~^\a) 

a. 

= N+e-** ^<5 CT , CT ,A CT , CT (q)e- 1 — , 

a, a' 

where we define 



and 



H = Hi + Hdi 



Hi 



(1) 



Pa> a(xm>\x m ) = ^™_f(ef - e Q ) (mV \a) (a\ ma) . 

a 

Here 5/ x ,, Xm +q v ) lii a periodic Kronig delta function, i.e. 
it is nonzero for x m > — x m + q y + s L x for any integer 
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s. With these definitions we can write the Hamiltonian 
matrix in the HF approximation in a compact form which 
is simple to diagonalize numerically 

(ma\n\m'a') HF = £ {( 7 A (q)U H (q) + U dis (q))I+ 

qeBZ 

|(A (q)I + A(q) • r)U F (q)} S^ Xm , +qy) e+^' 
-\gB-r, (2) 

where A Q (q) = Tr{A CTCT < (q)r Q }, and B specifies the ori- 
entation of the external magnetic field. The various ef- 
fective potentials which appear here are defined as 

Udis(q) = ^ e -ilq+G| 2 W£(q + G)e ife+G,)( 9H+ G y ) _ 

G 

(3) 

and 

U ^q) = -\ E e-i^e*^-^ ^(1 - V,o) , 

q' 

(5) 

wit hG= i^±n x ^n y ). 

For ue(t), we choose a white noise potential without 
spatial correlation, {{ve{*)ve ( r '))) = <J 2 8{r — r'). The 
density of states in the non-interacting limit has been 
calculated exactly Wegner for this distribution of the 
disorder potential,E-3 yielding a full width at half max- 
imum of approximately 1.06<r. In our calculations the 
parameter 7 = e 2 /ea specifies the relative strength of 
interactions and disorder broadening. This type of dis- 
order potential distribution is characterized by a single 
parameter a, which we use as our unit of energy. As we 
discuss later, our results are insensitive to correlations 
in the disorder potential on length scales smaller than I, 
but would change in some respects for disorder potentials 
which are smooth on the magnetic length scale. 

The HF equations are solved by an iterative approach 
which can create difficulties which must be addressed. 
The HF equations generally have many solutions that 
correspond to different, usually metastable, extrema of 
the HF energy functional. The challenge is to locate the 
true global minimum. In particular, the iteration pro- 
cess will not break any symmetries of the Hamiltonian 
which are not broken by the starting charge and spin- 
densities, even though the global minimum of the HF 
energy functional frequently does break at least some of 
these symmetries. To counter such problems, it is usu- 
ally a good idea to introduce small artificial terms in 



the Hamiltonian which break the continuous symmetries 
and help the iterative process to reach the lowest energy 
state. In this problem, the iterative process is also ham- 
pered by severe convergence problems at zero tempera- 
ture connected with the localization of HF quasiparticle 
wavefunctions and the long range nature of the Coulomb 
interactions. A small change in the energy of a partic- 
ular orbital may involve a substantial rearrangement of 
the charges. These problems can be mitigated by always 
working at a temperature which is comparable to the 
finite-size quasiparticle energy level spacing and which 
scales to zero as the system size increases. The Zeeman 
term in the Hamiltonian, for which we choose typical ex- 
perimental values, reduces the SU(2) spin symmetry to 
a U(l) symmetry. In order to break the continuous U(l) 
symmetry we introduce an artificial (but very small) local 
magnetic field at the center of our simulation cells which 
points in the x-direction. It is this space-dependent field, 
required on purely technical grounds, which has moti- 
vated developing the formalism in a manner which per- 
mits non-constant Zeeman coupling strengths. To ensure 
that this artificially field and the finite temperature do 
not affect the final solution, we lower the magnitude of 
these terms until no change is seen in local charge and 
spin densities, or in HF quasiparticle energies. 

The phase diagrams discussed in Sees. Ill and IV, 
were obtained by starting from the non-interacting case 
and incrementing the interaction strength 7, taking as 
the starting densities the self-consistent densities from 
the previous 7 value. There is, of course, some hysteresis 
involved in this process so we do a backwards sweep on 7 
once we have reached the maximum interaction strength 
for a given run. If they differ, we use the smaller of the 
values obtained in upward and downward sweeps for the 
energy per particle. The energy per particle is obtained 
using the expression 

f = ~£ U dis (q)A *(q) + 2- J2 U H (q)|A (q)| 2 

qSBZ qGBZ 

+U F (q) (|A n (q)| 2 + |A u (q)| 2 + |A u (q)| 2 + |A iT (q)| 2 ) 

~gB-P tot , 

where Ptot is the total global spin polarization. The local 
spin magnetization density, which we calculate as well, is 
given by 

(<%(») = % £ ^'(vNmK^CW^W . 

We define the local spin polarization as (Pi(r)) = 
2{Si{v))/{h{p{v))). Note that in this case (|P(r)|) does 
not have to be equal to 1 since the system is compressible 
except in the limit of very large 7 where (|P(r)|) — » 1 for 
allr. 

Our criteria for convergence is that 
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*a = f EE i A -'(q) - (q)i 2 < i >< io- 6 , 

* qGBZ aa' 

(6) 

where i stands for the ith iteration. We have performed 
calculations for several disorder realizations at different 
values of v for system sizes of = 16—32. The finite size 
effects come mainly from the effective exchange potett. 
tial Up(q) and have been studied in detail previously.E£l 
The main effect is on the interaction part of the en- 
ergy per particle and it is well understood and easily 
corrected. We believe that the physics of the phase tran- 
sitions observed in these calculations is not affected quali- 
tatively by finite-size effects. Our qualitative conclusions 
are based on persistent features which are obtained for 
several different disorder realizations. The values of 7 at 
which the various transitions and cross overs we discuss 
below take place, do not change by more than 5% for 
different realizations. The results we present here are for 
one particular disorder realization. 

III. RESULTS AT v = 1 

At v = 1 the disorder free (7 — > 00) 2DES has@ a 
S = N/2 ground state. The S z = S = N/2 member 
of this multiplet is a single Slater determinant and can 
therefore be obtained by solving Hartree-Fock equations 
self-consistently. It is only in recent years that samples 
which are sufficiently clean to reach, or at least nearly 
reach, complete spin polarization have been grown.tj The 
collective behavior producing such a ground state was 
not exhibited in earlier samples which had more disorder 
in the form of unintended impurities, interface disloca- 
tions, and, in modulation doped samples, the potential 
from remote ionized donors. Fig. |l| summarizes the HF 
theory results we have obtained for the dependence of 
the spin polarization on interaction strength. The cal- 
culations were performed for a realistic value of the Zee- 
man coupling strength, g = 0.015, and at a very small 
value, g — 0.0018. Extrapolating from these two val- 
ues to g = 0, allows us to identify parameter values for 
which spontaneous spin polarization occurs, i.e., values 
for which the ground state is ferromagnetic. We find 
that ferromagnetism occurs for 7 > 0.5 in the Hartree- 
Fock approximation; at smaller values of 7 the single- 
particle disorder term dominates and yields a spin-singlet 
ground state. Notice that the spin susceptibility, which 
may be estimated from the difference between the spin- 
polarizations at the two g values, is small in the singlet 
state, and becomes large as the phase transition to the 
ferromagnetic state is approached. For the specific finite- 
size disorder realization we have studied, complete spin 
polarization is reached at a finite value of 7 ~ 1.5. At 
larger values of 7, the system has a finite gap for charge 
excitations. We must be aware, however, that the HF 
approximation overestimates the tendency of the system 



to order so the interaction strength at both transition 
points should be taken as lower limits. In addition, any 
physically realistic disorder potential is likely to have rare 
strong disorder regions which prevent the fully polarized 
state from being reached. 

In our calculations, there is a wide region of interac- 
tion strengths 7 for which partially spin-polarized states 
occur. In this regime our HF ground states nearly al- 
ways have non-collinear magnetic order. We show local 
spin polarization and charge density profiles of typical 
partially polarized states in Figs. ^| and ||. The origin 
of the reduced spatially integrated spin polarization is 
partly due to variation of spin-orientation, but princi- 
pally due to a reduction in the average value of the of 
the magnitude of the local spin polarization. This point 
is illustrated in Fig. 1 (open circles) and may be inferred 
from Fig. || (a). The reduction in spin-polarization is 
due to the occurrence of doubly-occupied orbitals, i.e. 
to disorder induced charge fluctuations which cannot be 
accurately described in models which include only the 
spin degree of freedom. The charged excitations of the 
ground state in this regime are ungapped and involve 
population of localized quasiparticle states. We also re- 
mark that local density profiles at these relatively small 
7 values, illustrated in in Fig. || (a), follow the effective 
disorder potential smoothed by the form factor for lowest 
Landau level electrons. Rapid spatial variation compo- 
nents in the white noise model disorder potential have 
little effect on the electronic state. The relationship be- 
tween electron number density and the Pontryagan index 
density of the local spin orientation, valid for slow spirt, 
orientation variation and nearly constant charge densityEJ 
is not valid in this regime. Still, the collective nature 
of the 2DES manifests itself in the nonzero spin polar- 
ization density perpendicular to the Zeeman field. As 
interactions strengthen further, the local charge density 
smoothes out favoring the minimization of Coulomb en- 
ergy at a cost in disorder energy. (See Fig. 2 (b) and 3 
(b) for 7 « 1.5.) 

Experimentally, the effects of disorder can be seen most 
directly in the NMR spectral line shape obtained at the 
lowest possible temperatures where the_spin profile is 
frozen on the experimental time scale.E3E3 The NMR in- 
tensity spectrum in this regime is given by 

/(/, 7) oc J dr PA r(z)e-^( 27r/ - 2 ^^( z ^ 5 > ; ^» , (7) 

with a — 9.34ms _1 and K s ~ 25KHz. Here Pn(z) is the 
nuclear polarization density and p e (z) is the electron den- 
sity envelope function in the quantum well. The .evalua- 
tion of such spectra has been outlined elsewhere ;E3E3 here 
we simply show results for several interaction strengths 
in Fig. The parameters used in Eq. fj] are the same as 
the ones used in Ref. [l2]. 

Note that the quantity usually identified experimen- 
tally as the Knight shift, the location of the peak in the 
NMR spectrum in Fig. |[ does not match the global po- 
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larization. This Knight shift measurement always over- 
estimates the global polarization. In order to obtain the 
global polarization of the system from the measured spec- 
trum one has to extract the first moment of a normal- 
ized spectrum£3 One sees from the NMR spectrum at 
7 = 1.25 that disorder induced spin density variation 
leads to a broadening of the maximum peak and can even 
lead to secondary peaks at lower Knight shift frequencies. 
Note, however, that features corresponding to negative 
Knight shifts, corresponding to regions of reversed elec- 
tronic spins, are unlikely because the typical size of such 
regions is small and because they are also obscured by the 
finite width of the quantum wells which trap the 2DES. 
Our calculations demonstrate that care must be taken in 
interpreting low temperature NMR data in the quantum 
Hall regime. 

The partially polarized regime can also be studied ex- 
perimentally by measuring the transport activation gap. 
Provided that weak Zeeman coupling can be ignored, the 
extended quasiparticle states are expected to be precisely 
at the Fermi level when the 2DES is in a paramagnetic 
state. The Hall conductivity should jump from to 2e 2 /h 
at v = 1. In the ferromagnetic state, the majority-spin 
extended quasiparticle state will be below the Fermi level, 
the majority-spin extended state will be above the Fermi 
level and the Hall conductivity at v = 1 should be quan- 
tized at a xy = e 2 /h. The spontaneous splitting of the 
two extended state energies is experimentally accessible 
and should exhibit interesting power law critical behav- 
ior as the ferromagnetic state is entered. This transport 
gap should vary monotonically with the global spin polar- 
ization, although the precise relationship between these 
quantities is not trivial. 

It is interesting to note that the Skyrmion-anti- 
Skyrmiofl pairs predicted recently by Nederveen and 
Narazova do not appear in our calculations. We do not 
conclude that these objects cannot appear at v = 1; we 
would expect them, for example, if we choose a disor- 
der model with relatively large potential variations, but 
only on a length scale much larger than the Skyrmion 
size. In this case the NLer model considerations in Ref. 
H should be applicable. Our calculations demonstrate 
rather clearly however, that charge density variation at 
v = 1 does not necessarily, or even usually, require the 
existence of well defined Skyrmion quasiparticles. 

IV. RESULTS ATi//l 

In clean (large 7) samples where full polarization is 
observed at-p — 1 , the global polarization decays rapidly 
with 1 1 — v\ J3 It is generally accepted that this property is 
a unique signature which experimentally establishes the 
thermodynamic stability of Skyrmion collective quasipar- 
ticles. In the strong disorder limit, on the other hand, 
spontaneous spin polarization does not occur at any fill- 
ing factor near v = 1. 



The global polarization results for v 7^ 1 in Fig. [|, 
illustrate how the system interpolates between these two 
extrema. As the interaction strength 7 is increased from 
to 2, the behavior is similar to the v = 1 case. For 
strong disorder charge variation is dominant, and small 
spin polarizations occur primarily because many single 
particle orbitals are occupied by both up and down spin 
electrons. Charge variation is the dominant response to 
disorder, and it continues to play an important role at all 
interaction strengths. At sufficiently large 7, our finite 
size systems reach a state with the maximum spin po- 
larization allowed by the Pauli exclusion principle. This 
maximally polarized state is reached earlier than in the 
case at v = 1 (7 ~ 1.4—1.6) because, we believe, a 
larger number of charged quasiparticles are available to 
screen the random potential. At this point the system 
forms what we refer to as a conventional quasiparticle 
glass (CQG). The conventional Laughlin quasiparticles 
are initially localized in the deepest minima (or maxima 
for v < 1) of the effective disorder potential and as the in- 
teraction strength increases, or equivalently the depth of 
the disorder potential wells becomes smaller, the charged 
quasiparticles rearrange themselves locally into a quasi- 
triangular Wigner crystal pinned by the strongest of the 
disorder potential extrema. At larger 7 we observe a 
transition from a CQG to a Skyrmion glass. The loca- 
tion of this transition is marked by a reduction of the 
global polarization from its maximally polarized value. 
For a specific disorder realization the point of cross over 
from the CQG to the Skyrmion glass, as illustrated in 
Fig. [|, depends on filling factor and g. The dependence 
of the transition point on g in this regime can be ap- 
proximated by considering a simple model for a single 
Skyrmion trapped at a disorder potential extrema. We 
approximate its energy by 

E(K) = U(K - K a ) 2 + g*fi B BK + aAK , (8) 

where K is the number of spin flips per Skyrmion, a is 
the strength of the disorder potential and A is a phe- 
nomenological parameter. The first two terms determine, 
the optimal Skyrmion size in the absence of disorder .Ej 
The form for the third term reflects the property that 
Skyrmions with smaller K are smaller and will be able 
to concentrate more strongly close to the potential ex- 
trema. This simple model gives an estimate of the inter- 
action strength at which K > Skyrmions first become 
stable 



2K Q U/(e 2 /el) -g' w 

The parameters U and Kq can be estimatedEl for fill- 
ing factor v = 1.25 as U/{e 2 /ei) ~ 0.014 and K ~ 1. 
Using our numerical result for where the transition oc- 
curs at g = 0.0018, we estimate that A ~ 0.1. From 
this, one obtains an estimate of 7* ~ 7 for the cross over 
point from conventional quasiparticles to Skyrmions at 
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g = 0.015. This is in reasonable agreement with the ac- 
tual cross over point 7* ~ 10 (see Fig. |J the transition is 
out of scale in Fig. ||) given the simplicity of the model. 
These estimates of the maximum disorder strength at 
which Skyrmion physics is realized could be checked by 
performing NMR experiments in samples where electron 
density, and hence the interaction strength, is adjusted 
by the application of gate voltages. 

For a particular realization of the disorder potential, 
particle-hole symmetry is broken in a finite system and 
is recovered only in the limit of very large 7. The particle- 
hole symmetry relation for the global spin polarization is 
i.e. (1 - e)P z {u = 1 - e) = (1 + e)P x {v = 1 + e) where 
e < 1. At large 7 this relation is approximately satis- 
fied. Also in this limit, the Pontryagan relation between 
the local density profile and the local spin density£l be- 
comes accurate. In the clean limit, the Skyrmion system 
crystallizes in a square lattice for the filling factoxa con- 
sidered here. (The Skyrmion crystal is triangularis for v 
very close to 1.) The disordered Skyrmion glass state has 
very smooth fluctuations of the local spin density, com- 
pared to the CQG, although both lattices are pinned by 
the disorder potential. We remark that quantum fluctu- 
ations in Skyrmion positions are not accounted for in HF 
theory, and it is quite possible that even in this limit the 
ground state is a liquid rather than a crystalO As noted 
in Ref. [l(] it is possible that the broken U(l) symmetry 
of the Skyrmions orientation order predicted by Hartree- 
Fock does not survive quantum fluctuations. We show an 
example of the CQG in Fig. |] (a) and (a), and of the 
quasi Skyrmion lattice state in Fig. || (b) and [7] (b). Note, 
that we find, in agreement with Nederveen and Narazovjil 
a shrinking of the Skyrmion size as disorder broadening 
increases. This effect may help explain the appearance of 
a "tilted plateau" centered-around v = 1 in the Knight 
shift vs. filling factor dataJlH Rare highly disorder regions 
in the sample may localize and reduce the effective size 
of the few Skyrmions present at these filling factors. This 
would allow the bulk of the sample to be fully polarized 
at v 7^ 1 and give rise to a Knight shift equivalent to the 
one at v — 1 . The plateau is tilted because of the change 
in fully polarized density as pointed out in Ref. O. 

V. DISCUSSION 

We have used the Hartree-Fock approximation to study 
the competition between interactions and disorder near 
Landau level filling factor v — 1. At a qualitative level 
our results can be summarized by the schematic zero- 
temperature phase diagram shown in Fig. ^, which is 
drawn for the case of small but non-zero Zeeman cou- 
pling. Distinct ground states can be distinguished by dif- 
ferent values for the quantized Hall conductivity, a xy , by 
the presence or absence of spontaneous spin-polarization 
perpendicular to the direction of the Zeeman field, and 
by the presence or absence of a gap for spin-flip excita- 



tions. At small 7 (strong disorder), the electronic state is 
paramagnetic (denoted as PC in Fig. ||), there is no spin- 
polarization in the absence of Zeeman coupling, and the 
Hall conductivity is expected to jump from to 2e 2 /h as 
the filling factor v crosses the v = 1 line. For a small Zee- 
man coupling, there will be a small splitting between the 
majority-spin and minority-spin extended state energies 
and the zero-temperature Hall conductance should have 
a narrow intermediate e 2 /h plateau centered on v = 1. 
However, we do not expect that this plateau will be ob- 
servable at accessible temperatures, and have indicated 
this in Fig. || by using a thick line to mark the to 2e 2 /h 
phase boundary. At somewhat larger 7 there is a phase 
transition at zero Zeeman energy between paramagnetic 
and ferromagnetic states (denoted as FC in Fig. |J). I n 
our calculations this transition occurs at a larger value 
of 7 at v = 1 than away from v = 1. As 7 increases 
in the ferromagnetic state, we expect that the separa- 
tion between majority-spin and minority-spin extended 
state levels will increase rapidly so that the v = 1 integer 
quantum Hall plateau will broaden and become observ- 
able. At still larger 7, we find a transition to a state with 
the maximum spin polarization allowed by the Pauli ex- 
clusion principle. At v < 1, this is full spin-polarization. 
In these states, marked 'SC for spin-gap in Fig. |L the 
differential spin-susceptibility vanishes. For realistic dis- 
order models, it seems likely that in the thermodynamic 
limit there will always be rare high-disorder regions in the 
sample which prevent maximal spin-polarization from be- 
ing achieved. For this reason, the phase transition we 
find in our finite systems likely indicates a crossover from 
large to small differential spin-susceptibility in macro- 
scopic systems; we have therefore marked this transition 
by a dashed line. Finally at the largest values of 7 (weak- 
est disorder) the physics for v near 1 is dominated by 
Skyrmion quasiparticles which emerge from the v = 1 fer- 
romagnetic vacuum. In this regime, the system develops 
spontaneous spin-polarization in the plane perpendicular 
to the direction of the Zeeman field. In Fig. |] we have 
labelled this regime NCF for non-collinear ferromagnet. 

This phase diagram is intended to represent the fill- 
ing factor interval 0.85 < v < 1.15, over which fractional 
quantum Hall effects are not normally observed and it 
appears likely that Hartree-Fock approximation calcula- 
tions are able to represent interaction effects. Some of 
our findings may help explain the striking tilted plateau 
feature observed in the NMR spectrally near v = 1. Nev- 
ertheless, we have found rich structure in the crossover 
between non-interacting and disorder-free limits of the 
v = 1 quantum Hall effect which helps explain the diffi- 
culty experienced in attempting to construct a simple in- 
terpretation of low-temperature NMR spectra. Our cal- 
culations motivate experimental studies of v = 1 trans- 
port activation energy studies near the paramagnetic to 
ferromagnetic phase transition. 
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FIG. 1. Global polarization for a fixed disorder realization 
at v = 1 with g = 0.015 (filled circles) and g = 0.0018 
(squares). The average magnitude of the local polarization 
is shown (open circles) at g = 0.015. The transition to a fully 
polarized state occurs at 7 ~ 1.6 — 1.9 for all disorder real- 
izations studied. Examples of the local polarization density 
and local density profile at different values of 7 for a typical 
disorder realization are shown in Figs, fcq and H. 



7 



1 2 3 4 5 6 1 8 9 10 




FIG. 2. Local polarization density for v = 1 and g = 0.015 
at 7 = 0.75 (a) and at 7 = 1.25 (b). The z-component is 
represented by a shadow plot with black as — 1 and white as 
+1. The in-plane local polarization density is represented by 
a two-dimensional vector plot. The contour plot corresponds 
to the effective disorder potential used in the phase diagram 
shown in Fig. hi. 
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FIG. 3. Local density profile for v = 1 at 7 = 0.75 (a) 
and at 7 = 1.25 (b). In this shadow plot, black represents a 
local Landau level filling factor of two and white represents 
a local Landau level filling factor of zero. The contour plot 
corresponds to the effective disorder potential used also in the 
calculations of Fig. Minima in the electron density occur 
at maxima in the effective potential and vice versa. 
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FIG. 4. NMR spectrum for v = 1 and g = 0.015 at 7 = 0.75 
(solid line), at 7 = 1.25 (dashed thick line), and at 7 = 2.0 
(dashed-dotted line). The sample parameters correspond to 
ones the used in the experimental studies of Ref. [13] . 
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FIG. 5. Global polarization phase diagram for a fixed dis- 
order realization at v = 1.25 and g = 0.015 (solid line), at 
v — 1.25 and g = 0.0018 (dashed line), at v — 0.75 and 
g = 0.015 (dotted line), and at v = 0.75 and g = 0.0018 
(dotted-dashed line). The transition to a maximally polar- 
ized state (Laughlin quasiparticle glass) occurs at 7 » 1.5 — 2 
for all disorder realizations obtained. The transition from a 
Laughlin quasiparticle glass to a Skyrmion glass occurs at 
7 ~ 15 for experimentally relevant parameters. Local spin 
polarization density and local density at different values of 7 
for such disorder realization are shown in Figs. ^| and [?]. 




FIG. 6. Local spin-polarization density for v = 1.25 at 
7 = 10.0 (a) and at 7 = 40.0 (b). The z-component of the spin 
is represented by a shadow plot with black as —1 and white 
as +1. The in-plane local spin-polarization density is repre- 
sented by a two-dimensional vector plot. The contour plot 
shows the specific effective disorder potential which leads to 
these results and those illustrated in Fig. ^. Notice that spon- 
taneous in-plane spin-polarization appears only at the larger 
7 value. 
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FIG. 7. Local density for v = 1.25 at 7 = 10.0 (a) and at 
7 = 40.0 (b). The density is represented by a shadow plot 
where black correspond to local filling factor v = 2 and white 
corresponds to local filling factor v = 0. The contour plot 
shows the effective disorder potential used in this calculation 
and in Fig. ^. Notice that the density variation is smoother at 
larger 7 when Skyrmion quasiparticles, rather than Laughlin 
quasiparticles, are localized near potential minima. 
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FIG. 8. Phase diagram discussed in Sec. V. Here PC in- 
dicates the compressible paramagnetic phase, FC the par- 
tially polarized compressible ferromagnetic phase, FSG the 
spin gapped ferromagnetic phase, and NCF the non-collinear 
ferromagnetic phase. We emphasize that this phase diagram 
is qualitative in nature and transition points vs. 7 should be 
taken as upper limits. 
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